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Abstract 

Magnetotunnelling spectroscopy of resonant tunnelling structures provides information on the 
nature of the two-dimensional electron gas in the well. We describe a model based on nonequilib- 
rium Green's functions that allows for a comprehensive study of the density of states, tunnelling 
currents and current spin polarization. The investigated effects include the electron-phonon in- 
teraction, interface roughness scattering, Zeeman effect and the Rashba spin-orbit interaction. A 
qualitative agreement with experimental data is found regarding the satellite peaks. The spin po- 
larization is predicted to be larger than ten percent for magnetic fields above 2 Tesla and having a 
structure even at the satellite peaks. The Rashba effect is confirmed to be observable as a beating 
pattern in the density of states but found to be too small to affect the tunnelling current. 
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I. INTRODUCTION 



Resonant tunnelling was first observed in scattering of electrons by noble gases where it is 
known as the Ramsauer effect^. After the invention of quantum mechanics, the appearance 
of the scattering minimum at a particular electron energy was explained by coupling with 
a quasi-bound state of the gas atom. Following the pioneering work of Esaki and Tsu^ in 
the seventies, the same effect was found in resonant tunnelling structures^(RTSs) where the 
role of electrostatic potential barriers is played by thin semiconductor layers. 

RTS is the basic nanoelectronic device that exploits the quantum nature of electrons. As 
the electron transport through a RTS is a combination of coherent tunnelling and phase- 
breaking collisions, RTSs are the ideal testing ground for quantum transport theories. In 
the envisioned high-performance electronic applications, one aims to maximize the main 
current peak where coherent transport dominates while minimizing the satellite peaks that 
occur due to scattering. Resonances appear in scattering as well because it is enhanced 
for the localized quasi-bound states, the most pronounced feature being the longitudinal 
optical (LO) phonon satellite peak-. While leading to performance degradation and power 
dissipation, the scattering mechanisms in RTSs are interesting in themselves for allowing 
the study of various aspects of the twodimensional electron gas (2DEG) physics via cur- 
rent spectroscopy. When a perpendicular magnetic field is applied, the 2DEG density of 
states (DOS) collapses into discrete Landau levels (LLs) thus pronouncing the resonant fea- 
tures in scattering. For this reason, magnetotunnelling experiments in RTSs^~— have been 
established as an important tool in investigating the 2DEG physics. 

The electron-phonon interaction within the 2DEG of III-V RTSs leads to a formation 
of magnetopolarons^ manifested as a characteristic anticrossing of LLs which has been 
observed in magnetotunnelling experiments^-^. Another important phenomenon in InAs- 
based RTSs is the spin-orbit interaction (SOI). Due to potential applications in semicon- 
ductor spintronics^, the zero-magnetic field spin splitting in a 2DEG has been receiving 
considerable attention since the field effect spin transistor was proposed by Datta and Das^. 
At smaller magnetic fields SOI is known to significantly affect the 2DEG LLs leading to a 
beating pattern in the DOS. This causes a beating pattern in magnetoresistivity which is 
often used to characterize SOI in a 2DEG^~— . 

Here we report on a theoretical study of magnetotunnelling with LO phonon and interface 
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roughness (IR) scattering taking a particular care of the carrier spin. The spin splitting is 
described as a combination of SOI and the Zeeman effect, both of which are known to be 
pronounced in InAs quantum wells. A detailed description is given of a nonequilibrium 
Green's function (NEGF) model^— with scattering treated within the self-consistent first 
Born approximation (FBA). The LO phonon and IR correlation functions are approximated 
with delta functions, the suitability of which is discussed in Appendices |A] and El The 
model is applied to a symmetrical InAs-GaAs double barrier RTS and shown to yield I-V 
curves that qualitatively match the experimental data. An estimate of the spin polarization 
is found, and the role of SOI in spin-dependent transport through RTSs^i - — discussed. 

II. MODEL DESCRIPTION 

In the first part we discuss the electronic states of the RTS in absence of scattering. 
Eigenstates of the lateral and spin degrees of freedom found here are used in the second part 
to formulate the transport model with scattering. 

A. Effective Hamiltonian and the Landau levels 

The RTS is assumed to be grown along the [100] direction denoted as the z-axis. It 
comprises of: (1) a highly doped InAs emitter contact, (2) the 3nm wide left barrier formed 
by a GaAs layer, (3) the 6nm InAs well, (4) the 3nm wide right barrier made of GaAs and 
(5) a highly doped InAs collector contact. We adopt a one-band envelope function approxi- 
mation (EFA)^ in which the electronic states of a heterostructure are determined by their 
effective mass m* and the effective potential V(z) comprising the conduction band edge 
variation and the externally applied voltage. A generalized multiband EFA allows for many 
parameters and is applicable for modelling a wide range of semiconductor heterostructures^. 
The one-band approximation used here is often adopted for single-valley conduction band 
heterostructure states. This is because up to few hundred milielectron volts above the 
conduction band minimum, the states are of a dominantly s-like character so both energy 
dispersion and matrix elements can be sufficiently well estimated using effective parame- 
ters. The use of more complicated models such as the multiband EFA or, especially with 
NEGF calculations, multiband empirical tight-binding models^.^ is justified when a more 
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FIG. 1. The symmetric InAs-GaAs RTS under bias. The barrier and well widths are 3 and 6nm, 
respectively. The conduction band edge with slope due to the external bias is shown with the full 
(red) line. The shading represents the LDOS in logarithmic scale. The pale series of peaks in the 
emitter (left) and collector (right) region represent Van Hove singularities at LL subband bottoms 
while the intense peaks in the well are the quasi-bound LLs. Up to around 700meV only the lowest 
z-subband is present, the second z-subband states appearing as broadened double-node peaks. 

quantitative model^ is aimed for or even mandatory if effects like inter-valley scattering or 
inter-band tunnelling are studied. 

The conduction band effective potential V(z) of the investigated RTS is depicted by 
the thick solid line in Fig. ([TJ. The externally applied collector to emitter voltage U ce is 
assumed to be positive so that electrons have a tendency to flow towards the collector whose 
band edge is taken as a reference energy. The EFA Hamiltonian of the biased RTS with an 
external magnetic field applied parallel to the growth direction reads 

7T 2 1 OtR 

H = ^—^ + V(z) + -g*/j B Ba z + -^-{^ x a y - ir y a x ). (1) 

The first term is the kinetic energy, 7r = p + e A being the kinetic momentum of the electron 
and m* its conduction band effective mass. The third term represents the Zeeman effect 
determined by the effective ^-factor, g*, and (T X) y iZ are the Pauli spin matrices^. 

The fourth term is the Rashba spin-orbit interaction (RSOI) Hamiltonian that accounts 
for SOI due to the heterostructure inversion asymmetry-^. The Dresselhaus SOI (DSOI) 
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caused by the bulk inversion asymmetry is neglected. This is believed to hold in InAs 



quantum wells, considering the theoretical arguments in Ref. 
in Refs. 
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and 
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and experimental findings 
42 and 



43 (for InGaAs- 



(for InAs-GaSb quantum wells) and Refs. 
InAlAs quantum wells). The RSOI in a heterostructure is a combined effect of the externally 
applied potential and the spatial variation of the band edge. In asymmetric heterostructures 
the latter lifts the spin degeneracy even in absence of an external electric field E ext . For 
symmetric structures the Rashba constant or can be taken to be proportional to E ext 

a R = CE cxU C = C E + C h (2) 

where the direct contribution of the electric field is evaluated as** 

c _ eh 2 (2E g + A S o)A so 
E 2m* E g (E g + A so )(3E g + 2A so )' 1 J 

The interface contribution C\ depends on details of the heterostructure potential. We have 
found that in the range of external fields of interest (around the main and satellite current 
peaks) taking into account only the first term gives Q!r in the range 1 — 2 « 10~ n eVm which 
is close to values for a InAlAs-InGaAs quantum well reported that vary in the range 
2 — 5 ~ 10 eVm. For this reason, unless it is varied as a parameter, «r is evaluated from 
Eqs. (ED and (Ej) with £ g (InAs) = 356meV, A so (InAs) = 380meV while C\ is set to zero. 

Thus, our model has three main effective parameters, m*, g* and «r. The energy scale 
over which these parameters change considerably is set by the bandgap. In InAs the effective 
mass increases approximately linearly from the zone center value m* = 0.023m to around 
twice that at E — 0.5eV. The variation of the effective g-factor predicted by the k • 7r model 
is 



g*(E)=g*(0) 



E s (E g + A so ) 



A 



so 



E„ + E E B 



^so 



E 



with the zone center value g*(0) given by 



9*(0) 



m 



2A 



so 



(4) 



(5) 



m*(0) J 2 A so + 3£ g ' 

g*{E) calculated from the above equations with InAs parameters is a monotonously increas- 
ing function starting at g*(0) = —16.9 and going to around —4 at 0.5eV. The effective 
g-factor is also known to depend on the electron-electron interaction, thus a range of values 



has been reported in InAs quantum wells**. In Ref. |45| \g*\ of a 2DEG in a 4nm wide InAs 



quantum well has been found to be around 6. The situation with a R is less clear since it 
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FIG. 2. A scheme of the E ns levels and their relation with the ordinary (or = 0) LLs for negative 
g. (a) Shows the arrangement of the levels at high magnetic fieds. (b) is the dispersion E ns (B) for 
small and moderate magnetic fields. For clarity, (b) is drawn with exaggerated q;r. 

vanishes in bulk while a value in a heterostructure is highly dependent on both the con- 
stituents and electrostatic effects. It will be seen below that the mechanisms involved in 
electron transport through the investigated RTS have an energy scale of few tens of milielec- 
tron volts, determined by the chemical potential of the contacts and the InAs LO phonon 
energy. Consequently, keeping in mind the range of meaningful values, it makes sense to 
fix the effective parameters which is done by setting m*(InAs) = 0.023mo, g* = —5 and 
calculate as explained above. 

In absence of RSOI («r = 0), the eigenkets describing the lateral p = (x,y) and spin 
degrees of freedom are the spin split LLs \na) with energies 



En 



eB 

n = 0,l,..., A B = g*fi B B, u c = — , (6) 

m* 



where the plus is to be taken for spin up, a =t, and the minus for spin down along the z 
axis, cr —\.. Adopting the Landau gauge, the corresponding wavefunctions are 



[p | \kn t) = (pi \kn i) 



kL 



(pt \kni) = (pi |Jfcnt> =0. (7) 



In this equation k and L x denote the quasi-momentum and system length along the x axis, 
while L = y^H/eB is the magnetic length. The quantization of the lateral degrees of freedom 
is manifested in the local density of states (LDOS) shown in Fig. [T]where it leads to a series 
of Van Hove singularities (corresponding to the bottom of each LL subband) in the emitter 
and collector and to a series of quasi-discrete levels in the well. 

For nonzero RSOI the LLs given by Eqs. and Q are mixed into eigenstates labelled 
as | kns) where the allowed values of the pair of quantum numbers n and s are^— 

ns = 0+,l-,l+,2-,2+,.... (8) 

The corresponding energies are 



E ns =nhw c + s^ (^ hUJc 2 AB ^j +nR\ R=V2a R /L. (9) 
The | kns) states can be expressed in terms of \kna) as 

\kns) = Al\k(n - 1) t) + At\kn |), (10) 



with coefficients A a given by 



VnR i i{E n+ - E, 



n-ltJ 



y/nB? + {E n+ - E^)*' n+ y/nB? + (E n+ - E^)*' 

Al = , , At = - & (11) 

^nR^iE^-E^ ^JnR? + (E ni - E n J) 2 

In evaluating Eq. fflOl) for for ns = 0+, \k(n — 1) t) is to be taken as the zero ket. The 
arrangement of the E na and E ns energy levels for a negative effective ^-factor is depicted 
in Fig. [2J By analyzing Eqs. ($\) and ffTUl) . several observations can be made regarding the 
effect of RSOI on LLs : (1) the |0+) ket is identical (apart from the irrelevant phase factor) 
to \0 1) and they have the same energy, E 0+ = E ^, as depicted by the dots connecting these 
two levels in Fig. [2j (2) the two closely spaced spin-split levels E n ^ and E n ± become E( n +x)- 
and E n+ , respectively; (3) the Zeeman spin splitting, A B , becomes AE n = E n+ — E( n +i)- 
which increases with n and is greater than A# for all n. 



B. The quantum transport model 



Initially, magnetotunnelling in RTS has been studied^— within the frame of the 



Landauer-Biittiker formalism neglecting the phase 
ing ad hoc schemes to describe it^. In Refs. [52 



Dreaking scattering altogether or us- 



5J the Matsubara technique suitable 



for a rigorous treatment of many-body effects at finite temperatures^ is used to evaluate 
the RTS transmission probability in the presence of LO phonon scattering. The combined 



electrostatic charging and LO phonon scattering effects have been studied in Ref. |56| by 
evaluating a conveniently devised equation of motion for the electron field operators. While 
the degree of complexity and ability to account for various physical phenomena of the mod- 
els used in these studies vary, their common difficulty is a lack of a consistent description 
of the non-equilibrium state of the biased RTS. It is for this reason that the NEGF model 
that simultaneously accounts for the quantum mechanical details and the non-equilibrium 
statistics is the principal method in studying the RTS. Existing NEGF studies of RTS 
magnetotunnelling are scarce and include the study of scattering on LO phonons^ and the 
effect of shallow impurity states^. 

The quantities used in the NEGF method are the Green's functions G 1 and self-energies 
S 7 both of which can be of three kinds: lesser (7 =<), greater (7 =>) and retarded 
(7 = R). In explaining the model we need to jump from the mixed, \kns)\z), to the real- 
space, I per) 1 2) = |rcr), representation and vice versa, so to improve the paper readability we 
start by giving a step-by-step explanation of how the equation for determining G R is arrived 
at. 

The retarded Green's function G R is defined in the time domain as^ 

G R lS2 (k l n 1 z u k 2 n 2 z 2 ]t 1 M) = - l -6{t 1 -t 2 ) ({$ nfamzn h), ^\k 2 n 2 z 2 ] t 2 )}) . (12) 

Here ti% are the time coordinates, 9(t) is the Heaviside step function, (...) is the ensemble 
average and {...} the fermion anticommutator. ty s (knz, t) and ty' s (knz, t) are the destruction 
and creation operators for the \kns) state at coordinate z and time t. The spin coordinate 
(s = ±) is arbitrarily written as a subscript. As the problem under consideration is station- 
ary, all the quantities depend only on the time difference t\ — 1 2 and we switch to the energy 
domain, e.g. 

G^Jhrnz^fam^E) = j d(h -^e^'^G^Jhn^h^z^h ~h). (13) 

Since the Hamiltonian ([T]) allows for decoupling the lateral and the z degrees of freedom while 
I kns) are eigenstates, G R is diagonal in kns in absence of scattering (any perturbation to H). 
The onset of a scattering mechanism generally mixes the \kns) states so off-diagonal terms 
may appear. However, depending on the scattering model details, the coupling between 
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any given \kin\Si) and all the other \kns) states induced by the perturbation can often 
be properly described by a diagonal term of the self-energy E B . All the Green's functions 
and self-energies considered here are diagonal in k and n but not in s, so the full notation 
from Eq. f lT5]) will be abbreviated by G^ lS2 (knziz 2 ; E) . This (quasi-) decoupling of the kn 
coordinates will be clarified later after the scatterer correlation functions are introduced. 
The retarded Green's function is found by fixing kn and inverting the matrix 

Gf lS2 (kn Zl z 2 ; E) = (E- H SlS2 {kn Zl z 2 ) - Zf 1S2 {kn Zl z 2 ; E))' 1 , (14) 

where the z degree of freedom is represented by a tight-binding model with the nearest- 
neighbour hopping energy^ t = h 2 /2m*a 2 and care taken to ensure Hermiticity around the 
heterointerfaces where the effective mass changes abruptly^. It is equivalent to replacing 
the derivatives along z by finite-differences defined on a set of equidistant points with a 
being the point-to-point distance. Therefore, if the z axis interval enclosing the RTS is 
represented by N mesh points, Eq. (]14p is a 2N x 2N matrix equation, the doubling being 
due to spin. 

The self-energies in (|T4"|) are a sum of four terms that are calculated self-consistently 

^ = E^ + E£ + E? R + £Z . (15) 
The emitter (E^) and collector (E^) self-energies describe the interaction of electrons within 



the RTS with the exterior of the device, as explained in Ref. [26j. The interface roughness 
(E7 R ) and LO phonon (E^ R ) self-energies are assumed zero in the first iteration. 

The lesser (A =<) and greater (A =>) Green's functions are found from the kinetic 
equation 

dz 3 / dz4GK S3 (knziz 3 ;E)Y,* 3S4 (knz 3 z^E)Gf 4S2 (knz4Z 2 ;E), 

(16) 

where G A is the advanced Green's function (in the energy domain, it is the complex conjugate 
of G R ). In subsequent iterations FBA is used to evaluate the interface roughness (<fi = IR) 
and phonon (0 = LO) lesser and greater self-energies 

^>ir 2 ; E) = JL J ^G^(nr 2 ; E - E^D^r,; E% (17) 

where D^(tit 2 ; E) denotes the scatterer correlation function assumed to be independent on 
the electron spin (we neglect the scattering via SOI). For IR scattering the lesser (-D < ) and 



greater (D > ) correlation functions are equal 

D m (nr 2 ; E) = 2 7 rh5(E)5(r 1 - r 2 ) £ E&(zO<f(* - (18) 

where the summation is done over the heterointerfaces located at z — Z\. The LO-phonon 
correlation functions are taken to be 

L>g ( ri , r 2 ; u) = 2nh6(E T £ L o)0( r i - r 2 ), (19) 

where E^q = 29meV is the LO phonon energy in InAs. The delta-like form of correla- 
tion functions given by Eqs. 018p and (fTUj) is used because it considerably simplifies the 
calculations^ 1 ^ by decoupling the equations for states with different kn coordinates. To 
clarify the approximation being made, note that the IR correlation function used in the 
literature is usually a Gaussian function^ with the deviation set by the roughness lateral 
length A (typically few nanometers), while the Frohlich interaction with LO phonons gives a 
dependence proportional to^ l/|ri — r 2 |. Roughly, the approximation consists in assuming 
that the self-energy at a given point in space is determined solely by the electronic prop- 
erties at that same point while the more realistic models allow for a contribution from all 
the nearby points. A quantitative discussion of scatterer correlation functions aimed at de- 
termining meaningful values for scattering strengths Ujr and {7lo is given in Appendices lAl 
and El 

The mixed representation of the lesser and greater self-energies is then found as 
EjL iaa (*™> *i> ^ w) = S( Zl - z 2 ) UUz»i - *i) E G ^(Pi^u P^ 2 ; E)A°* sl A" nS2 , (20) 

I cr 

for IR scattering and 

^ Qsi sM *i, *2| u>) = 6( Zl - z 2 )Ul £ G%{ Px z 1 , Pl z 2 - E =f E hQ )AZA° nS2 , (21) 

(j 

for LO phonons. The real-space representation of Green's functions in the above equations 
is given by 

GZAfhzr, PxZxM = —^ J2A° Sl A'Z 2 Gl S2 {n,z 1 ,z 2 ;u), a=t,|- (22) 

ns 1 s 2 

The retarded self-energy is related to the lesser and greater self-energies by 

E R (E) = \ [£>(£) - E<(E)} + ^-VV J dE ^ >{E 2Z^ (E '\ (23) 
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where VV denotes the Cauchy principal value. The left term is the anti-Hermitian part that 
causes the broadening of energy levels. The right term is Hermitian and corresponds to the 
scattering-induced energy shift (the energy correction found by time-independent perturba- 
tion theory^). In studying the electron transport through a RTS, the energy shift (of the 
order of a milielectron volt) due to various interactions is not very significant especially since 
the exact quasi-bound state energies are hard to estimate anyway. On the other hand, the 
level broadening leads to qualitative differences and its exact value is crucial in estimating 
the degree of spin polarization. Therefore, in order to avoid the difficulties in the numerical 
evaluation of the Cauchy principal value, the Hermitian part is neglected, as often done in 
NEGF calculations^. 

Once the values of T, R are found, the calculations are repeated to obtain a self-consistent 



solution, as in Ref. 
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III. RESULTS AND DISCUSSION 



Magnetotunnelling through a RTS is governed by the coupling of the emitter and collector 
electron baths via the discrete states in the well, the coupling being a combined effect of 
coherent tunnelling through the barriers and incoherent transport assisted by scattering on 
the interface roughness and LO phonons. Consequently, the starting point of our discussion 
is the (quasi) discrete LL spectrum in the well, which is described by the spin-polarized 
DOS in the RTS. After clarifying the structure of this spectrum, we move on to study 
the properties of the vertical transport and discuss the features in the RTS I-V curves. 
Throughout this section low temperature (T = 4.2K) with equal emitter and collector 
electrochemical potentials (/ig = A*c — 20meV) are assumed. 



A. Properties of the RTS DOS 

In absence of RSOI, the dispersion of LLs, E ncr (B), is represented by the usual Landau 
fan diagram. The presence of RSOI leads to a nonlinear dispersion E ns (B) given by Eq. 
OH]) and plotted for first few LLs in Fig. [2j By investigating Eq. ([9]) it is seen that the 
arrangement of E QS levels shown schematically in Fig. [2] is reached at B large enough that 
fo^c 3> y/nR. The large spacing between the doublet E n+ and E( n _^_ and any other level 
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implies a simple DOS structure consisting of a series of double peaks (separated by the 
spin-splitting energy AE n = E n+ — £( n _i)_) arranged at a distance hue apart. However, a 
rich structure is seen in DOS at smaller magnetic fields. 

Considering that fiwc is proportional to B and that R is proportional to \fB, at B — > 
all the E n _ levels lie below all the E n+ levels (Eq + being the lowest amongst them). With 
increasing B, each E n _ level approaches the corresponding £'( n+1 ) + level crossing all the 
E n+ levels below it. This produces a characteristic beating pattern in the DOS for any fixed 
B^-. While the exact energy spacing AE between the two lowest nodes is easily calculated 
from the E ns (B) dispersion, to understand its order of magnitude and dependence on B and 
the RSOI parameter we note that the crossing occurs when the energy difference hwc 
between two successive LLs is compensated by ^JnR. Therefore n ~ (hwc/R) 2 and 

AE ~ nhw c = (fkucf/R 2 , (24) 

where the tilde denotes the order of magnitude. As an example, for an = 0.936 x 10~ 11 eVm, 
B = 0.5T (giving huj c ~ 2.39meV), Eq. ([24]) predicts AE ~ 99.7meV. 

The spin-polarized DOS in the RTS is obtained by integrating the spin-polarized LDOS 
over the length of the RTS 

DOS(£) CT = / dzLDOS a (E, z), (25) 

J Z E 

with 

LDOS CT (£, z) = (G* (pz, pz- E) - Gt(pz, pz; E)) , (26) 

and the total DOS is a sum of its spin polarized components, DOS(-E) = DOS^(-E) + 
DOS;(£). 

There are two conditions that need to be satisfied for observing the crossing of the few 
lowest LLs: (1) R ~ Tvjjq and (2) hue > T, where T denotes the LL broadening. The latter 
condition ensures that B is large enough for a discrete spectrum to appear. In a RTS, T is 
nonzero even for noninteracting electrons due to the finite probability of electron escaping 
from the well via tunnelling through barriers. For the structure shown in Fig. (TJ T = 0.8meV 
in the ballistic limit (scattering neglected). Therefore, the LLs become clearly resolved for 
fields around 0.5T (hue ~ 2.5meV). At external voltages around the main current peak 
(U ce = 270mV), from Eq. ([3]) we find «r = 10~ 11 eVm, giving R ps 0.4meV, so the crossing 
would not be observable. For purposes of illustration, Fig. [3] shows the spin-polarized DOS 
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FIG. 3. Spin-polarized DOS of the RTS biased with U cc = 270mV showing the dispersion E ns (B). 
To observe the LLs in this range of magnetic fields, the LLs must be very narrow. Here the 
ballistic limit is taken while the value of ckr is exaggerated to five times the value predicted by Eq. 
[3] (approximately 4.6 x 10~ 11 eVm). The plotted quantity is log [1 + DO S a {E) / Z>2DEg]- 

of the RTS for external bias U cc = 270mV (main peak of the current) calculated in the 
ballistic limit and with «r set to five times the value given by Eq. ([3]) so that the crossing 
of LLs is clearly seen. The 2DEG spin-polarized DOS at zero magnetic field is given by 
-D2DEG = m*/2nh 2 . The energy is always measured relative to the collector band edge, so 
that the emitter band edge is U e = eU ce . If the RTS had a DOS as shown in Fig. [3] it would 
mean that the RSOI could be investigated by studying the structure of the RTS I-V curves. 
However, considering that the realistic value of V is few meV and that the effective mass in 
an InAs quantum well is likely to be higher than 0.023mo, we conclude that it is not realistic 
to expect that the RSOI effects would affect the shape of the I-V curves. 

When scattering is included with strengths £/j 2 R and Ul Q as found in Appendices lAl and IB"l 
the broadening for magnetic fields around IT becomes T « 3.5meV. It causes the overlap- 
ping of individual LLs for smaller magnetic fields and, consequently, makes the observation 
of the nonlinear features of the E ns dispersion impossible. However, RSOI can be experi- 
mentally observed via the beating pattern in Shubnikov-de Haas (SdH) oscillations^— . It 
is a consequence of the beating pattern seen in the DOS. This is illustrated in Fig. HI The 
beating pattern survives even with significant scattering effects. The distance between the 

13 




DOS/D„ „ 

2DEG 

FIG. 4. DOS of RTS biased with U ce = 270mV and B = 0.5T. The red and black curves were 
calculated for a moderate IR scattering intensity (A = lnm) . The difference between the two is the 
onset of LO phonon induced broadening at energies more than Slo above the emitter conduction 
band bottom. In SdH measurements this effect is not observed since the LO phonon emission at 
the Fermi level is blocked by the Pauli exclusion principle. The blue curve is the case with higher 
IR scattering, showing that the beating pattern exists even when the LL broadening is larger than 
fuj c . 

nodes is found to be around lOOmeV, in good agreement with the estimate given by Eq. 

(El. 

The effect of LO-phonon scattering becomes apparent if the DOS is observed on a larger 
scale. Figure |5] shows the (a) ballistic and (b) LO-phonon scattering DOS of the biased 
RTS. Compared to Fig. [5] (a), additional energy levels are seen in (b). These are the phonon 
replicas of the LLs displaced from them by multiples of £lo- The problem of electron-LO 
phonon interaction in a 2DEG, being of fundamental significance, has received a considerable 
attention^^S. Polaron effects have been studied using cyclotron resonance^ in a quantum 
well, or by investigating the I-V curve features in magnetotunnelling in perpendicular-*^ and 
in-plane^ magnetic fields. 

The DOS shown in Fig. [5] (b) demonstrates that the current model captures some of the 
main features of the electron (magnetopolaron) DOS in a (quasi)2DEG: (1) appearance of the 
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(a) ballistic (b) LO phonon 




2468 2468 



B [T] B [T] 

FIG. 5. The evolution of DOS with B for B > 1. The RTS is biased with U cc = 270mV. (a) 
Ballistic limit where the LL broadening is entirely due to coherent coupling with emitter and 
collector states, (b) DOS when the LO phonon scattering is included with U^q = AUq . The red 
line shows the emmiter Fermi level. The Pauli blocking of the first replica of the lowest LL level 
below 4T (when it is below the red line) is not complete because the system is not in equilibrium 
so the lowest LL is only partially occupied. The plotted quantity is log [1 + DOS(E) / D2deg]- 

phonon replicas, (2) the anticrossing of LLs and their replicas at magnetophonon resonances 
and (3) a considerable broadening of LLs only around the magnetophonon resonance. Effects 
related to the real part of the self-energy, see Eq. ( I2"3"j) . such as the polaron shift and effective 
mass renormalization, are not included in this model. The LO phonon scattering strength 
C/lo determines the height of the DOS corresponding to phonon replicas. It also affects the 
anticrossing energy. A quantitative connection between the LO-phonon self-energy and U^q 
is not simple to establish, owing to the fact that the self-energy is found as a self-consistent 
solution, but a rough estimate is that the self-energy scales as the square root of Ul Q . The 
significance of the phonon replicas in the RTS is that they provide an additional channel for 
electrons tunnelling through the RTS. For example, the lowest spin-split replicas seen in Fig. 
E](b) correspond to the lowest lying |0+) and |1— ) LLs. Once their energies are aligned with 
electrons in the emitter having the same lateral coordinates, ns = 0+ and ns — 1— , resonant 
intra-LL (preserving the LL coordinates) tunnelling with LO phonon emission occurs. 

The effect of IR scattering on the DOS is simpler and consists of broadening all the 
states (both the LLs and their replicas). In contrast to the phonon scattering, it is known 
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B=5T 




D0S a( E ) /D 2DEG 

FIG. 6. The effect of IR on spin-polarized DOS, DOS^- is shown in full lines and DOS^ in dashed 
lines. The bias voltage is U cc = 270mV, while the magnetic field is B = 5T. 

to depend only on the DOS and not on the occupation of the levels 2 ^. In Appendix |A] we 
argue that if a delta model is used for the roughness correlation, the scattering strength at 
interface z\ should be taken as 

U m = V h 2 A 2 1R 7rA 2 . (27) 

In this approach only the product of Air and A is significant, so we fix Air = 0.3nm and 
consider the effect of changing A. Values of A = 5 — 7nm have been reported for AlGaAs 
quantum wells^ 1 ^ and yield a broadening of more than 15meV. Since the typical level 
broadening in a high quality RTS is expected to be several^ meV, smaller values of A have 
been used in calculations. In Fig. [6] a comparison is made of the spin-polarized DOS in 
the biased RTS at B = 5T for the ballistic (A = 0), A = lnm and A = 2nm cases. The 
IR scattering increases the width of the LLs to around 4meV for A = lnm and 8meV for 
A = 2nm. The LL profile is found to significantly deviate from the Lorentzian profile which 
can be explained by the discrete spectrum (a Lorentzian profile is obtained if the retarded 
self-energy is approximately constant around the given level which is a characteristic of a 
continuous spectra). 
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B. Terminal currents 



The paramagnetic current flowing through the RTS is obtained by evaluating the trace 



of the current operator I op (E). Following Ref. 
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, instead of I op , we evaluate 



KpW = ^ [X<(E)G>(E)-Z>(E)G< (E)] , p = E,C,LO (28) 

which has the same trace as I op (E). The symbol p denotes the system with which the current 
is exchanged^. For calculating the I-V curves, the emitter (E) and collector (C) currents 
are needed. The exchange of particles with the phonon bath is described by the p = LO 
term. Utilizing the decoupling of kn coordinates, K p (E) in the mixed representation reads 

K PsS1S2 (knzi, z 2 ; E) = / dz 3 [H< slS3 (knz 1 ,z 3 ;E)G> 3S2 (knz 3 ,z 2 ;E) 

-j:> jS1S3 (knz l ,z 3 ;E)G< S2 (knz 3 ,z 2 ;E)] . (29) 

The diagonal elements K P)Cra (pz, pz; E) in the real-space representation are found from a 
transformation given in Eq. (1221) . Finally, the spin and energy resolved current density is 
obtained by integrating over z 

i Pt(T (E) = J dzK Ptaa (pz,pz;E). (30) 

The total currents are obtained by integrating over all the energies 



dEi p>a (E). (31) 



Charge and spin conservation implies z'e,ct = ic,a, so if only the total currents are required, 
either of the two can be calculated. However, the emitter and collector currents have a 
different energy distribution when dissipation (due to inelastic scattering) is present which, 
as shown below, provides information on what happens with electrons on their way through 
the RTS. 

We discuss the details of electron transport through RTS on the example of B = 5T. 
While the choice of a particular value of B is not critical, this one has been taken based 
on several criteria. Firstly, it is large enough {FujJq = 23.9meV) that only the lowest LL 
subband of the contacts is populated, which makes the discussion simpler. This is because 
the electrochemical potential /i of both contacts at B = is assumed to be p = 20meV, while 
at B — 5T it is found to be around p, = 18.2meV (also see the Fermi level variation in Fig. 
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FIG. 7. The energy resolved spin-up collector current ic^{E) evolution with U ce , low voltage 
range. To emphasize the low value peaks, the map shows log [1 + 100 x iq <[(E)/I max \, where 
/ mffl = 114.6A/cm 2 meV is the highest value of ic +(E) in the investigated range. The inset shows 
the total current ic + where the horizontal axis is U cc in units of mV and the vertical axis is given 
in units of ^4/cm 2 . The enumerated peaks are explained in the text. 



[5]). The dependence of /i on B corresponds to a three dimensional system (the contacts are 
also assumed to be exposed to the magnetic field). Secondly, the LL spacing is larger than 
the level broadening, so the features in the I-V curves can be observed clearly. Thirdly, the 
lowest magnetophonon resonance (denned by E^q = hue) is close so we can discuss how it 
is manifested in the RTS currents. Finally, the Zeeman energy Ab = — 1.45meV at this field 
determines the spin-splitting (of the lowest LLs) while RSOI is practically negligible. The 
consequence is that the LLs are polarized approximately either spin-up (\kn+) « \kn \)) or 
spin-down (\k(n + l)—) m \kn f)) so the transport of the two spin polarizations is practically 
decoupled. We will consider the spin-up current which is dominated by the \kn—) states, the 
admixing with \kn+) states being negligible. For this reason, in explanations we mention 
only the \kn— ) states, even though the mixing is taken into account in numerical calculations. 

The transport of electrons across the RTS is studied by considering the energy resolved 
collector spin-up current ic^(E) shown in Fig. [7J Red lines show the position of the emitter 
band edge U e (lower), LL subband bottom U c + hw^jl (middle) and the top of the emitter 

18 



Fermi sea, U c + fi. The zero values between U e and U e + huc/2 are due to that fact that the 
emitter DOS is zero in this energy range. We do not consider the emitter current ie^(E) 
separately, as it offers little information. This is because we know by default that all the 
electrons that leave the emitter have energies between U e + huc/2 and U e + \i. The total 
current is shown in the inset. For U cc < 240mV the currents are practically zero because 
the quasi-bound states are above the emitter Fermi sea. 

The peaks in iq^(E) labelled by 1, 2 and 3 for U ce between 260 and 290 form the main 
peak of the current. In this range of U ce the emitter Fermi sea is aligned with the lowest lying 
| AT— } quasi-bound state, see Fig. 1 denotes the electrons that preserved their energy 
while going through the device. The majority of this current represents the electrons that 
have coherently tunneled through the RTS while a small part of them have been elastically 
scattered into any of the LL subbands in the collector. 2 is exactly E^o below 1, implying 
that it represents electrons that have entered the well, emitted one LO phonon and got 
scattered in the collector. Similarly, 3 is the current component of resonant transmission 
with emitting two LO phonons. 

The second peak of the total current seen in the inset of Fig. [7] occurs when U ce is 
increased enough so that the \k2— ) quasi-bound level and the phonon replica of \kl— ), 
see Fig. (b), are aligned with the emitter Fermi sea. The peak labelled by 4 represents 
resonant transmission mediated by IR scattering of } levels into the \k2— ) levels in the 
well. It is seen as a shoulder of the peak in the inset. In agreement with Fig. [5] (b), 4 occurs 
for slightly lower U cc than 5, which represents the resonant transmission via the phonon 
replica of \kl— ). With increasing B, 4 and 5 are found to move in opposite direction, 4 
moving towards higher U cc . The anticrossing of the LLs is seen as a disappearance of 4 at 
the magnetophonon resonance (around 6T) and reappearance at higher U ce , so that 4 and 5 
never occur at the same value of U ce . Resonant transmission with two LO phonon emission 
is again observed as 6. The resonant enhancement of this process is a consequence of the 
increased population in the ) level into which electrons get by emitting the first LO 
phonon. Some of the electrons then pass coherently to the collector (5) while the remaining 
emit an additional phonon. 

The same situation but for U ce ranging from 360 to 480mV is shown in Fig. [HJ The most 
apparent difference compared to the lower voltage range is the appearance of higher order 
(up to five) LO phonon emission processes. This is because with increasing U ce more LLs 
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FIG. 8. The energy resolved spin-up collector current iq ^(E) evolution with U ce , the high voltage 
range. The map shows log [1 + 100 x ic \{E) / I max \, with I ma x = 3.2A/cm 2 meV. The inset is the 
total current Iq + where the horizontal axis is U ce in units of mV and the vertical axis is given in 
units of A /cm 2 . The enumerated peaks are explained in the text. 



in the well are available to scatter into by emitting an LO phonon. The peak labelled by 7 
occurs when \k3— ) is alligned with the emitter Fermi sea. At slightly higher U ce appear 8 
and 9 corresponding to the first replica of \k2— ) and second replica of \kl— ), respectively. In 
the high voltage range an important difference between the elastic and inelastic currents is 
seen: the elastic current appears almost as a structureless stripe in the U e + tvjJc/2 to U c + fi 
energy range. The peak labelled as 10 corresponding to IR scattering via the |fc4— ) state is 
hardly visible. This is probably due to the gradual shifting of the LLs in the well towards 
the collector, which decreases their interaction with emitter states via the IR of the emitter 
barrier. On the contrast, the inelastic current comprises many peaks, each representing 
the alignment of a phonon replica with the emitter Fermi sea. 11 corresponds to the third 
phonon replica of the \kl— ) state. When considering the entire range of voltages, it is seen 
that the peaks 5, 9 and 11, corresponding to the first, second and third replica of \kl— ) are 
the most pronounced in the inelastic current. It seems counterintuitive that 9 and 11 should 
be higher than some of the peaks with less phonons emitted, because the two-phonon or 
three-phonon emission processes are weaker that the one-phonon emission process. This is 
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FIG. 9. (a) I-V curves for with B changed as a parameter. In (b) the I-V curves for U cc > 300mV 
are shown magnified by a factor of 10. The curves in (a) and (b) are vertically displaced by 
500A/cm 2 and 50A/cm 2 , respectively. The red lines represent the spin-up, while the black lines 
are the spin-down currents. 

also seen in Fig. |5](b) where the second replica of |M— ) clearly has a lower DOS than the first 
replica. While the results indicate that second and third order processes make a difference 
in calculations, it should be noted that the reason that the two- and three-phonon peaks 
in 9 and 11 are more intense then the lower order peaks is the sequential emission of one 
phonon at a time. This is possible, even though the magnetophonon resonance E^o = 
is not exactly achieved, due to the significant level broadening caused by IR scattering. 

The effect of varying the magnetic field on the total currents is shown in Fig. [91 There is 



very little difference between the B = 1' 
Fig. [9] (see, for example, Fig. 4 in Ref. 



and B = 0T curves, the latter not being shown in 



30l ). The is because the broadening of quasibound 



states is larger than frujc = 4.8meV. For B = 2T, the emitter chemical potential is around 
20meV so two LL subbands are occupied, leading to the structure in the main peak (the 
onset of the coherent \k2— ) current is around U ce = 265mV), and the satellite peaks - the 
two peaks in the 310 — 340mV range are the first replicas of \kl— ) and \k2— ). At B = 3T, 
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FIG. 10. Spin polarization P for various B. The curves are vertically displaced by 0.2. 

\i = 21.6meV, while 3/2 x hue = 21.5meV, so the second emitter LL subband is almost 
cut off. The coherent current from the \k2— ) emitter subband participates in the current 
but cannot be observed in the I-V curve. From B = 4T and above, only the lowest LL 
emitter subband is occupied. The peaks seen in the I-V curves occur when the alignment 
between either a higher LL or a phonon replica with the emitter Fermi sea is achieved. The 
evolution of the peaks towards higher values of U ec reflects the evolution of the DOS with 
B, as discussed on the example of the B = 5T case. 
The spin polarization P of the current, defined as 



with the magnetic field varied as a parameter is plotted in Fig. [TOJ The highest value 
of P is always reached at the main peak but significant values occur at satellite peaks as 
well. To explain this, we argue that the value of P is determined by two factors: (1) the 
magnetization of the emitter Fermi sea which becomes significant once only the lowest LL 
subband is populated and (2) the spin splitting of quasi-bound LLs in the well. Obviously, 
the magnetization is also determined by the spin splitting (but that of the emitter states). 
However, the distinction is important because (1) and (2) represent two mechanisms leading 



P 



2C,t - ^CU 



(32) 
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to nonzero P: the former influences the electron supply while the latter determines how well 
will the RTS itself filter the electrons based on their spin. P is higher at the main peak 
because most of the current is coherent and the quasi-bound LLs involved are narrower. 
However, the fact that the structure in P reflects so clearly the scattering processes means 
that (2) is significant even at the satellite peaks. 

IV. SUMMARY 

A comprehensive model for quantum transport in RTSs in a perpendicular magnetic 
field has been described. Effects like the formation of magnetopolarons and IR induced LL 
broadening are accurately reflected in the evolution of the I-V curves with changing the 
magnetic field. Using effective parameters that realistically reflect the InAs bandstructure, 
a beating pattern in the DOS has been obtained even with large IR scattering induced 
broadening, which is in agreement with experimental findings that the RSOI constant can 
be estimated from SdH measurements. It has been found that the I-V curves are not 
sensitive to RSOI, and that the spin polarization is dominated by the Zeeman effect. The 
spin polarization of the current is found to be significant for magnetic fields of few Tesla 
even at satellite peaks, suggesting that its measurement might yield additional information 
in magnetospectroscopic studies. 
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Appendix A: Correlation function for IR scattering 

The quantity used to describe interface roughness is the interface fluctuation ApR(r), 
i.e. the length by which the interface is displaced from the average value. By definition, 
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the average (over r) value of ApR(r) is zero. Typical values of Air are of the order of 



one monolayer, i.e. few angstroms (see e.g. Ref. l68l ). therefore the ^-dependence of the 
scattering potential caused by the roughness, Vni(r), can be taken as a delta function and 
Air considered to be a function of only the lateral coordinates Ai R (p) 

V m (r)=5(z-z l )V h A lR (p), (Al) 

where Vt, represents the conduction band discontinuity at the heterointerface. 

If the interaction between electrons and a scatterer is described by the Hamiltonian 

H' = J dr^{r;t)V{r;t)^{r;t), (A2) 

where V(r; t) can be a static external potential (in which case t is redundant) or involve 
dynamic degrees of freedom as in case of phonons, the scatterer correlation functions are 
defined as^ 

D > (n, r 2 ; t h t 2 ) = (F(n; ti)F(r 2 ; t 2 )) , D<(n, r 2 ; t h t 2 ) = (V(r 2 ; t 2 )V{n; h)) . (A3) 

The brackets (...) imply the statistical average which is the configurational average in case 
of the IR scattering or the ensemble average in case of phonons. 

For IR scattering the interaction potential Vi R (r) is static so the lesser and greater cor- 
relation functions are equal. After transforming to the energy domain, we find 

D m (r 1 ,r 2 ;E) = 2vh5{E)5{z x - z I )6(z 1 - z 2 )V*(A m ( Pl )A m (p 2 )}. (A4) 

The quantity (Ai R (p 1 )Ai R (p 2 )) is called the IR autocorrelation function in the literature 
and usually assumed to be of the Gaussian forrn^i^ 

(A rR (p 1 )A IR (p 2 )) av = Ae~^) 2 /A 2 ; ( A5) 
where A is referred to as the roughness height and A as the roughness lateral length. The 



values reported in Ref. [68| are A = 3 — 5A and A = 50 — 70A In problems where A is small 
compared to other relevant lengths, the Gaussian may be approximated by a delta function, 

so 

D(r ir2 ; E) = 2ivh5(E)5(z 1 - Zl )5{ Zl - z 2 )5{ Pl - p 2 )V*A 2 nA 2 . (A6) 

When a magnetic field is perpendicular to the interface the magnetic length L = ^JhjeB 
is the relevant quantity. As it becomes smaller than lOnm for B > 6.2T, we find that the 
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delta approximation cannot be considered as strictly accurate in the investigated range of 
magnetic fields (up to 10T) but that it can still be expected to give qualitatively correct 
results. 

Appendix B: LO phonon correlation function 

For polar coupling to bulk LO phonons (Frolich interaction), the term V(r; t) in Eq. flA2j) 
is given by^ 

^) = X> q e-*M qW , M^l^(-±- r -L), (Bl) 

where Aq(t) = a q (t) + a q (t) is the phonon operator while a q (t) and a q (t) are the phonon 
destruction and creation operators, respectively. Assuming that the phonons are in a ther- 
modynamic equilibrium and neglecting the LO phonon dispersion, = E^q, we find 

(A^A^h)) = 5 qi ,_ q2 {e-'^ofr-ftO/fi [n B (E LO ) + 1] + e^^^^V^o)} • (B2) 

At low temperatures such that k^T <C E^q, the Bose-Einstein occupation factor is very 
small, n-B(E^o) <C 1, so it can be set to zero (only spontaneous phonon emission is con- 
sidered). Upon switching to the energy domain and summing up over q, the LO phonon 
correlation functions are found as 

D^ (r u r 2 ;uj) = 5(uTUlo) \ — r , C LO = ( _J_ - -}—) . (B3) 

|ri-r 2 | Aeq \e{oo) e{0) J 

In the main text, the correlation functions are approximated by 

£>! LO (ri, r 2 ; w) = 2ir§(u t w LO )^Lo 5 ( r i - r 2)- (B4) 

In contrast to the case of IR scattering where the scatterer correlation is a short range 
Gaussian function, approximating the long range 1/r correlation by a delta function may 
appear bizarre. The first issue to be resolved when the 1/r correlation is used is finding 
the relevant length scale. In extended systems, such as the bulk, the scale is set by the 
screening length. As the problem of electron-phonon interaction in a RTS is analogous to 
the problem of electron-phonon interaction in a quantum well, we may consider the latter, 
leaving screening aside, and conclude that for our purposes the scale is set by the quantum 
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well width L z and the magnetic length L. But since L z is the extension of the entire 
investigated system, approximating D^ Q by Df LO is obviously flawed in a fundamental way. 
However, we will now show that the scattering strength L^o can be chosen so that Df LO 
gives scattering rates for electronic states in a quantum well that are close to scattering rates 
obtained using D^ Q . 

We consider an infinitely deep quantum well of width L z in a perpendicular magnetic 
field, k, n and the space coordinates are the same as in the main text. As RSOI and the 
Zeeman effect are irrelevant, we neglect them. The eigenkets are \kna) where a denotes the 
quantum number for motion along the z axis. As the aim is only to compare the two models, 
it is sufficient to consider FBA without self-consistency. The Green's functions below are, 
therefore, the non-interacting Green's functions with only diagonal elements in kna, labeled 
by G^(na; E). The self-energies for the delta model are evaluated straightforwardly as 

^Ilo^i™! !' hn 2 a 2 ; E) = 5 klk2 5 nin2 ^j^ G^(n 3 a 3 ; E =f E lo ) / dz \(z\a 2 )\ 2 \(z\a 3 )\ 2 . 

(B5) 

The self-energies for D^ Q are not diagonal in kn, but for the first order effect we need to 

26 



consider only the diagonal elements labelled by T^ Q (kna;E) 

^ (kna;E) = J2 G *( n ^ E T E LO )j^ j d^Ml\F(q z ,a x ,a)\ 2 \H{q h n u n)\ 2 , (B6) 

where 



niai 



F(g z ,«i,a)= I dz(a\z)e~ zqzZ (z\a 1 } 1 (B7) 

and 



|iJ(<?H,ni,7 



n! /£ 



2 \ ni—n 

2 = p-« 2 /2 " 



j ni — n 



ni > n, £ = q\\L. (B8) 



In these equations q is the phonon wavevector with qy and q z being its in-plane and per- 
pendicular components while is the generalized Laguerre polynomial. For more details 
on evalua ting th e electron-LO phonon matrix elements^ and scattering rates in a 2DEG 
i|. Assuming that only the lowest lying a = 1 subband in an infinitely deep 



see Refs. 
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quantum well is populated (as is the case with the RTS investigated in the main text) with 



{z\l) = y/2/L~ z sm(7rz/L z ), the self-energies given by Eqs. ( 1B5I) and (1B6|) are seen to be 
equal if 



Ulo = L2 ffl EL ° f^-T " -L) Q(nx,n 2 ,B), (B9) 

f |F(g z ,l,l)| 2 |^(gil,"-i,"-2)| 2 
Q(n 1 ,n 2 ,B) = J dq l W " 11 1 2,1 . (BIO) 

The value of Ul Q given by Eq. (1B9I) ensures that the scattering rate from LL rii to the 
lower lying n 2 due to LO phonon emission obtained assuming the delta correlation Df LO is 
the same as the one obtained from the exact D^ Q . The dependence of Ul Q on B and LLs 
involved in scattering is analyzed by noting that only n\^n 2 pairs which differ in energy by 
E^o are relevant. Therefore, we fix the lower LL index n 2 in Eq. (IB9j) and numerically 
integrate Q(ni, n 2 , B) to obtain the scattering strength as a function of B, while choosing 
the upper LL index rii so that rii = n 2 + [E lo /}tw c ], where the square brackets [u] denote 
the integer closest to u. The calculated L^o f° r L z = 6nm are shown in Fig. HH in units of 
Uq = lmeVnm/D2DEG- These results show that the scattering strengths that should be used 
to accurately describe the strength of various inter-LL transitions in the investigated range 
of magnetic fields are well within the order of magnitude with each other. Consequently, 
a Df LO correlation with a fixed (i.e. independent of ni,n 2 and B) value of U~l may be 
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expected to give predictions similar to those of a D^ Q model. 
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